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ABSTRACT 

We report on an investigation of the SBS 1520+530 gravitational lens system. We have used 
archival HST imaging, Keck spectroscopic data, and Keck adaptive-optics imaging to study the lensing 
galaxy and its environment. The AO imaging has allowed us to fix the lens galaxy properties with 
a high degree of accuracy when performing the lens modeling, and the data indicate that the lens 
has an elliptical morphology and perhaps a disk. The new spectroscopic data suggest that previous 
determinations of the lens redshift may be incorrect, and we report an updated, though inconclusive, 
value ziens = 0.761. We have also spectroscopically confirmed the existence of several galaxy groups 
at approximately the redshift of the lens system. We create new models of the lens system that 
explicitly account for the environment of the lens, and we also include improved constraints on the 
lensing galaxy from our adaptive-optics imaging. Lens models created with these new data can be 
well-fit with a steeper than isothermal mass slope (a — 2.29, where p oc r~") if Hq is fixed at 72 km 
s~^ Mpc""'^; isothermal models require Hq 50 km s~^ Mpc~^. The steepened profile may indicate 
that the lens is in a transient perturbed state caused by interactions with a nearby galaxy. 
Subject headings: galaxies: individual (SBS1520-I-530) 



1. INTRODUCTION 

The strong gravitational lens system SBS 1520-1-530 
(hereafter SBS1 52 0) was first investigated by 
IChavushvan et aTl (|1997f) . who found that the sys- 
tem consists of a pair of images of a broad absorption 
line quasar {zgrc = 1.855) separated by 1'.'6. The lensing 
galaxy was soon detected wit h adaptive optics (AO) 
imaging (jCrampton et al.|[l998D and was assumed to be 
at the redshift of one of two absorption line systems 
seen i n the spectra of the quasar images. iBurud et al.l 
(|2002f ) attempted to deconvolve the lens spectrum from 
the quasar spectra and found the redshift of the lens 
to be ziens = 0.717. This redshift is broadl y consistent 
with a photometric redshift determined by IPaure et al.l 
(|2002f ). Furthermore, the lens was found to lie along 
the line of sight to a photometrically identified cluster 
of galaxies that is expect ed to be at approximately the 
same redshift as the lens (jBurud et al.ll200 iFaure et"aT] 
[2002). 

Optical monitoring campaigns have led to the mea- 
surement of a time delay between the quasar images o f 
- 130 days (|Burud et al.ll2002t iKhamitov et al.ll2006f ). 
This time delay provides an additional constraint for 
determining the mass slope of the lens galaxy (e.g., 
iRusin et al.ll2003l ) and allows a value of the Hubble Con- 
stant to be determined for the system {Hq = 51 km s~^ 
Mpc~^assuming an isothermal mass profile: [Burud et al] 
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I2002f ). Note, however, that the mass slope and, thus, 
the value of Hq depend on the environment surround- 
ing th e lens system (e.g.. lDobke et al.ll2007l : [Auger et al.1 
l2007f ). An incorrect understanding of the mass distri- 
bution and environment of the lens might account for 
the anomalously low value of Hq obtai ned for SBS1520 
compared to other len s systems (e.g., iKoopmans et al 



2053| : I York eFaH 120051) and the WMAP (iSpergel et al 



20071) and Hubble Sva ce Telescope (HST) Key Project 



(jFreedman et al.l[2001h results. 

In this paper we present new Keck AO and archival 
HST imaging of the lens system that indicates the lens- 
ing galaxy may have a disk component. We also present 
a spectroscopic investigation of the lens environment and 
provide a new analysis of the lensing galaxy's spectrum 
which results in an updated lens redshift of zi^ns = 0.761. 
We discuss the implications of these new observational 
data on previous analyses performed with SBS 1520 and 
suggest that, in spite of some complexity, this lens pro- 
vides an interesting platform to investigate dark matter 
interactions in dense environments. 

2. IMAGING 

High resolution imaging of SBS1520 at 2.12 /im {Kp- 
band) was obtained on 2006 May 2 with the Near In- 
fraRed Camera 2 (NIRC2; Matthews et al. in prepa- 
ration) behind the adaptive-optics bench on the W. 
M. Keck-II Telescope. The NIRC2 narrow camera was 
ucd£li3g.i3d&>r the imaging, which provides a field of view of 
X 10" and a pixel-scale of 9.94 mas pixcl^^. 
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data were taken in nine 90-s exposures with a small dither 
between each exposure to facilitate good sky background 
subtraction during the reduction process. A sodium laser 
guide star was used to correct for the atmospheric tur- 
bulence, and a, r ~ 12 star, 15" from SBS1520, was used 
for Tip Tilt correction. 

The data were reduced within iraf^ using a double- 
pass reduction algorithm. The first pass is used to de- 
termine the positions of the objects on the sky, which 
are then masked out during the second pass where the 
usual sky subtraction with temporally adjacent frames is 
carried out. Bad pixel and cosmic ray masks are gener- 
ated during this process and these masks are then used 
to construct weight maps. The individual exposures are 
combined using th e weight maps and the drizzle package 
(|Fruchter fc Hoo9 12002) . which also corrects the data 
for the geometric distortions across the NIRC2 camera 
to provide precise relative astrometry. The point-source 
full-width at half maximum in the final drizzled image is 
54 milliarcseconds. 

The inner 8" of the AO imaging is shown in Figure [U 
clearly showing the two lensed quasar images, the lens- 
ing galaxy, and a nearby galaxy. The high resolution 
and sensitivity of the observations allow us to determine 
the morphological parameters of the lensing galaxy with 
a high degree of confidence. The A and B images of 
the quasar were fitted and subtracted using an empirical 
PSF determined from the star to the southeast of the lens 
(Figure [2K); the K-band image of the lens shows an early 
type morphology, perhaps with a disk. A best-fit galaxy 
model was then determined for the lensing galaxy using 
a non-linear fitting code to model the lens with a Sersic 
model convolved with the empirical PSF. The residuals 
from this model are shown in Figure [21d; the best-fit Ser- 
sic index is n = 2.7 and the effective radius is determined 
to be Te — 0'.'49. The ellipticity from the model is found 
to be 0.44 at a position angle of 156°. 

We have also obtained HST WFPC2 and NICMOS 
archival data for SBS1520 {HST GO Programs 8175, 
7495, and 7887; PI Falco). The imaging data consist 
of 2100 s with the F555W filter and 1600 s with F814W 
using WFPC2, and 2816 s with F160W and the NIC2 
camera. Th ese data were reduced using the multidriz- 
zle^ package (|Koekemoer et al.ll2002l ) , and galaxy proper- 
ties w ere measured using SExtractor (jBertin fc Arnouti 
Il996l) . The relative positions of the sources used in our 
analysis have been measured from the NIRC2 imaging 
if the objects are in the NIRC2 field of view; otherwise, 
the F814W imaging has been used. The lens system 
was placed at the center of the PC chip for the WFPC2 
data and in the center of the NICMOS field, providing 
complete coverage out to 16" for the WFPC2 fields (see 
Figure [31) and 12" for the NICMOS field. The object 
~ 1"8 to the Northeast of the lens is not resolved in con- 
ventional ground-based imaging but is revealed to be a 
galaxy in AO and HST imaging, with a FWH M approxi- 
mately 30% larger than the PSF. Following Burud et al.l 
([2b02) ■ we will denote this galaxy as Galaxy M. 

^ IRAF (Image Reduction and Analysis Facility) is distributed by 
the National Optical Astronomy Observatories, which are operated 
by AURA, Inc., under cooperative agreement with the National 
Science Foundation. 

^ multidrizzle is a product of the Space Telescope Science Insti- 
tute, which is operated by AURA for NASA. 



3. SPECTROSCOPY 

We have obtained spectroscopic observations of 
SBS1520 usin g the Echellett e Spectrograph and Im- 
ager (ESI: .Sheinis et al.ll2002l ) on the Keck II telescope 
and the Low-R esolution Imaging Spectrograph (LRIS; 
lOke et anil995f ) on the Keck I telescope. Our primary 
lens spectroscopy was obtained with ESI on 2000 July 3 
in conditions with thin cirrus and approximately 0"8 see- 
ing. Four 1800 s exposures were observed with a 0'.'75 slit. 
The ESI slit was oriented to run through both quasar 
images and through the lensing galaxy. Two LRIS slit- 
mask observations were made on 2004 June 15 in non- 
photometric conditions with approximately 1"5 seeing. 
The purpose of the LRIS observations was to measure 
redshifts for both the lensing galaxy and nearby galax- 
ies. The 831/9200 grating, with dispersion 0.9 A pixel"^, 
was used on the red arm of the spectrograph and the 
blue arm employed the 600/4000 grism with dispersion 
0.63 A pixel- ^ Four exposures of 1800 s were obtained 
for each mask, though one exposure of the second mask 
was discarded due to poor transparency. Both of the slit- 
masks contained a slit placed over the lens system in the 
same orientation as the ESI slit, providing seven good 
LRIS exposures of 1800 s for the lens system. 

3.1. Data Reduction 

The LRIS data were reduced using a custom auto- 
mated pipeline written in Python (Auger, in prepara- 
tion). The pipeline removes the instrumental bias and 
overscan regions of the LRIS CCDs, corrects for ampli- 
fier gain offsets, and flat fields the data from the red 
arm of the spectrograph. The program then automati- 
cally determines a wavelength solution from arclamp ex- 
posures and skylines, corrects for telluric absorption fea- 
tures with an airmass-scaled model of the typical absorp- 
tion at Keck, perfo rms an 'optimal' background subtrac- 
tion (jKelsonI [2003( 1 . and resamples and coadds the two- 
dimensional spectra to a constant wavelength-scale grid. 
The pipeline also searches each slit for spectral traces or 
emission lines and performs an aperture- weighted extrac- 
tion. The extracted spectra were cross-correlated with 
galactic and stellar template spectra to determine red- 
shifts for the science targets. All of these redshifts were 
then verified by hand; several were discarded due to poor 
quality spectra while others were corrected by manually 
identifying features. 

The ESI data were reduced with a modified version 
of the Python scripts used to reduce the LRIS data. 
A star was observed at several positions along the slit 
and the traces of the star were used to straighten the 
echelle dispersion orders. A wavelength solution was 
determined from HgNe, CuAr, and Xe arclamp obser- 
vations. The science data and an ob servation o f the 
spectrophotometric standard star IIZ44 (IOkd[1990[ ) were 
then straightened and the night sky background was sub- 
tracted. The spectrum of IIZ44 was then extracted and 
an atmospheric absorption and instrumental response 
model was created. The SBS1520 quasar spectra were 
then traced and extracted. The galaxy spectrum was 
extracted by defining a Gaussian aperture with a full- 
width at half maximum of one pixel centered slightly 
off of the position of the lens galaxy to minimize the 
shot noise from the quasar light. The extracted spec- 
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Fig. 1. — Keck NIRC2 adaptive optics K-band image of the SBS1520 lens system. The image is centered on the lonsing galaxy, and two 
field stars are indicated. The galaxy to the north east of the lens is denoted Galaxy M. The lensed quasar components are also indicated. 
The resolution of the image is O'.'OS. 



tra from the four exposures were then coadded and a 
smoothed version of the galaxy spectrum was created 
by convolving the galaxy spectrum with a Gaussian ker- 
nel of width (T = 100 kms^^. We also extracted the 
galaxy spectrum by modeling t he quasar ima ges and lens 
gal axy in a method simi l ar to iBurud et all (12002') (also 
see ICourbin et all 120001 : lEigenbrod et al.l l2007t ). Two 
stars (shown in Figure [T]) were observed on the ESI and 
LRIS slits along with the lens system to help define the 
PSF, but wc found that our observations did not have 
sufficient sensitivity to adequately model the galaxy and 
quasar spectra without substantial artifacts. 

3.2. SBS1520 Quasar and Lens Spectra 

Previous observations of SBS1520 have revealed sev- 
eral line-of-sight absorption systems present in the quasar 
spectra at z Ki 0.7 17 and z k, 0.815 (Chavushvan ct al. 
[19971: IBurud et al.ll2002^ . iBurud et all (|2002f) found evi- 
dence for Call K & H absorption features associated with 
the z — 0.717 system in the two quasar spectra and also 
find these lines in a spectrum of the lens galaxy obtained 
by deconvolving the lens spectrum from the quasar spec- 
tra. The absence of Call features associated with the 
z = 0.815 system led those authors to assign the lens 
galaxy the redshift z = 0.717. Our data confirm the ab- 
sorption systems seen at z = 0.717 and z — 0.815, but 
the quasar spectra also reveal a Mg and Fe absorption 
system at z = 0.760 and a CIV absorption system at 
z = 1.728 (Figure [3]). Furthermore, we have also found 
evidence for a galaxy at z — 0.761 that might be the 



lensing galaxy, with the z = 0.717 galaxy being a local 
perturber. The three lower redshift absorption systems 
are potentially associated with the lensing galaxy, and we 
evaluate the likelihood that the lens is at each of these 
redshifts. 

3.2.1. z = 0.815 Candidate 

Previous studies did not find evidence for a stellar pop- 
ulation associated w ith the z = 0.815 absorption system 
(|Burud et al.ll2002t l and we are also unable to find any 
broad absorption features that might be associated with 
the stellar population of a massive lensing galaxy; there 
were no features associated with the z = 0.815 system 
in our extracted galaxy spectrum. The quasar spectra 
reveal a Mg and Fe absorption system that also con- 
tains narrow Call H & K absorption lines, though these 
Call features are not associated with a stellar popula- 
tion. Furthermore, we note that the absorption lines are 
stronger in the spectrum of quasar image A than quasar 
image B, which is contrary to expectation because the 
galaxy is much closer to image B (see Figure [1]). A tab- 
ulation of the absorption features at z = 0.815 is found 
in Table [T] The absence of a stellar component and the 
stronger features being present in the spectrum of im- 
age A make it highly unlikely that the primary lensing 
galaxy is associated with this absorption system. Strong 
Mgll absorption systems have been ob served at high im- 
pact parameter (jC hurchill et al.l l2005l and there are sev- 
eral galaxies within 16" (approximately 80 h~^ kpc at 
2 = 0.815) of the lens that could be responsible for the 
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Fig. 2. — Keck AO K-band image of SBS1520 with (a) the quasar images subtracted and (b) the quasar images and the Sersic-model of 
the lensing galaxy subtracted. Circles indicate the centers of the subtracted components; note that a faint core is noticeable in the residuals 
of the galaxy. 
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Fig. 3. — WFPC2 F814W image of the SBS1520 field. All galaxies within 16" of the lens are indicated with circles and have been 
explicitly included in the lens modeling. The bright star 15" west of the lens was used for the adaptive optics imaging. 



quasar absorption lines. 

3.2.2. z = 0.717 Candidate 

The SBS1520 quasar spectra also show very clear evi- 
dence of Call absorption at a redshift z = 0.717 (Table 
[U associated with a system of Mg and Fe absorption. 
The structure of this absorption system's Mgll lines at 
2796A and 2803A is complex and contains features that 
span a velocity range of approximately 550 km s~^; this 
large spread in the velocities of the substructure compo- 
nents might be indicative of infall in a forming galaxy. 
The spectrum from the B image of the quasar does not 
show as much substructure as the A component, though 
the velocity spread is essentially the same as the A com- 
ponent (Figure E]). The absorption lines are much more 
prominent in the spectrum of quasar image A compared 
to quasar image B; this is counter to what might be ex- 
pected if the absorption system is associated with the 
lens galaxy because the lensing galaxy is closer to image 
B. 

The spectrum of the lens galaxy, extracted from an 
aperture between the quasar spectra, displays the broad 
Call absorp t ion fe atures at z — 0.717 identified by 
iBurud et al.l (|2002| ). The hues are weak (see Figure 
but are at the same redshift as the absorption syst em. 
The presence of these lines led iBurud et alJ ()2002D to 
ascribe the redshift of this system to the lensing galaxy. 
However, evidence against this identification includes the 
stronger absorption features in the A component com- 
pared to the B component, the large amount of substruc- 
ture seen in the A component indicative of a forming 



galaxy, and the presence of a stronger set of stellar ab- 
sorption features at a redshift z — 0.761. 

3.2.3. z — 0.761 Candidate 

Another pair of broad Call lines are found in the 
galaxy spectrum at 6986.64A and 6926.06A (z = 0.7606, 
Figure O. The Call H feature from this system has a 
line strength nearly two times greater than the Call H 
line from the z = 0.717 system. The line at 6986A might 
al so be seen in the s pectrum of quasar B in Figure 7 
of IBurud et al.l (|2002[ ). though we do not find evidence 
for the Call lines in our spectra of either of the quasar 
images. There is an Fe and Mg absorption system dis- 
covered in our ESI and LRIS quasar spectra (Table [3]) at 
the redshift za = 0.7603. The absorption features in the 
spectrum of quasar image B have greater rest equivalent 
widths than the corresponding features in the spectrum 
of image A (Table [3]), suggesting that the absorption sys- 
tem is closer to image B than image A. There is also a 
velocity offset of vb-a ~ 200kms~^ between the absorp- 
tion seen in the A and B images of the quasar; this sug- 
gests that there may be two galaxies associated with the 
cluster at z = 0.76 (see Section [373t responsible for the 
quasar absorption spectra or there may be a rotational 
component associated with the possible disk seen in the 
AO imaging (Figure [2^). 

3.3. Lens Environment Spectra 

The two LRIS slitmask observations of the field of 
SBS1520 have yielded redshift identifications for 52 
galaxies. We treat the lens galaxy as being at a redshift 
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Fig. 4. — There are four absorption systems seen in the SBS1520 quasar spectra; a portion of the extracted spectrum from the A 
component is shown above. Prominent features from the z = 0.717 (dash-dotted), z = 0.760 (dotted), and z = 0.815 (dashed) systems are 
indicated; most of these features are Fe and Mg hnes. The spectrum has been divided by the continuum and several of the Mgll lines are 
seen to be saturated. 



TABLE 1 

QSO Absorption at . 



: 0.8154 



Component A 



Component B 



Line 


REW (A) 


Wavelength (A) 


REW (A) 


Wavelength (A) 


FeII2344 


0.58 


4255.448 


0.34 


4255.589 


FeII2374 


0.29 


4310.385 


0.15 


4310.473 


FeII2382 


1.20 


4325.458 


0.75 


4325.630 


FeII2586 


0.52 


4695.580 


0.30 


4695.703 


FeII2600 


0.87 


4720.156 


0.69 


4720.290 


MgII2796 


1.68 


5076.165 


1.35 


5076.384 


MgII2803 


1.42 


5089.242 


1.09 


5089.421 


MgI2852 


0.41 


5179.019 


0.25 


5179.208 


CaII3934 


0.12 


7142.672 






CaII3969 


0.14 


7205.966 







z = 0.761 a nd use the objectiv e group finding algorithm 
described in lAuger et al.l ()2007D to find groups in the field 
of SBS1520. To summarize that algorithm, we define a 
maximum rest-frame velocity difference between galaxies 
in a candidate group; galaxies with velocity differences 
larger than this are placed into a new candidate group. 
We then define a mean group center and redshift for each 
group candidate, determine a velocity dispersion for the 
group using the gapper method, and apply an iterative 
clipping algorithm to remove outliers. This procedure 
leads to the identification of three spectroscopically con- 
firmed groups in the SBS1520 field, each associated with 
one of the absorption systems described in Section 13.21 



In Table [3] we describe the group properties, including 
the mean position, the number of galaxies in the group, 
the offset of the group centroid from the lensing galaxy, 
the redshift, velocity dispersion, and virial radius of the 
group, and the convergence (k) produced by the group 
at the location of the lens assuming a singular isothermal 
sphere (SIS) mass profile for the group. Table [S] contains 
information for the galaxies that populate each group. 

The spatial distribution of these systems is shown in 
Figure [7l Note that our data have a paucity of spec- 
tra in the immediate vicinity of the lens due to several 
factors. First, both of the slitmasks that we used in 
the observations contained a slit at least 10" long placed 




Fig. 5. — The Mgll doublet at 2796A (dashed hnes) and 2803A (solid lines, offset by 0.1 for clarity) for the z = 0.717 absorption system. 
The lines from the B component have been shifted upward by one unit for clarity and the zero point for the velocity has been taken to be 
the center of the deepest feature in the spectrum of the A image. 



TABLE 2 
QSO Absorption at z = 0.7165 



Component A Component B 

Line REW (A) Wavelength (A) REW (A) Wavelength (A) 



FeII2344 


0.80 


4023.84 


0.25 


4021.76 


FeII2374 


0.42 


4075.75 


0.07 


4073.63 


FeII2382 


1.24 


4090.04 


0.59 


4088.17 


MnII2576 


0.18 


4423.17 






FeII2586 


0.70 


4440.07 


0.30 


4437.73 


MnII2594 


0.23 


4453.47 






FeII2600 


1.38 


4463.29 


0.65 


4461.24 


MgII2796 


2.58 


4799.89 


1.49 


4797.73 


MgII2803 


2.13 


4812.28 


0.94 


4810.11 


MgI2852 


0.69 


4897.08 


0.13 


4894.67 


CaII3934 


0.45 


6753.93 


0.08 


6750.38 


CaII3969 


0.25 


6813.67 







over the lensing galaxy and oriented along the axis join- 
ing the two quasars, precluding the placement of slits 
across immediately neighboring galaxies. Additionally, 
the masks were approximately centered on the lens and 
therefore the LRIS slitmask bar obscures galaxies near 
the lens system. Finally, both masks placed a slit over a 
very bright star near the lens to help with point spread 
function modeling instead of placing the slit on neighbor 
galaxies. Therefore, the redshifts of galaxies very near to 
the lens were not sampled and this may systematically 
increase our estimate of the group offsets from the lens 
center (and thereby decrease the estimates of k). 



4. DISCUSSION 

4.1. The SBS1520 Lens Redshift 

The appearance of two sets of stellar absorption lines 
in the quasar spectra indicates that the lens must be 
modeled with at least two galaxy components (e.g., 
iBurud et al.ll200"^ . However, it is not clear at which 
redshift the lensing galaxy is located or which of the 
neighboring galaxies could be assigned the redshift of 
the second absorption system. The presence of stronger 
absorption in the B component compared to the A com- 
ponent for the z — 0.761 system, the relative strengths 
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Fig. 6.— ESI spectrum of the lensing galaxy of SBS1520 between 66OOA and 7200A. Call K & H from the lens galaxy at z = 0.761 
are marked with dashed lines; features from an intervening absorption system at 2 = 0.717 are marked with dotted lines. The spectrum, 
smoothed with a Gaussian kernel of width a = 100 kms"^, is overplotted. 



TABLE 3 
QSO Absorption at z = 0.7603 



Component A Component B 



Line 


REW (A) 


Wavelength (A) 


REW (A) 


Wavelength (A) 


FeII2344 


0.07 


4126.442 


0.08 


4129.185 


FeII2382 


0.18 


4194.293 


0.25 


4196.938 


FeII2586 


0.08 


4553.237 


0.07 


4556.341 


FeII2600 


0.21 


4577.002 


0.14 


4580.098 


MgII2796 


0.38 


4922.277 


0.55 


4924.809 


MgII2803 


0.31 


4934.881 


0.42 


4938.293 



of the stellar Call features between the z — 0.717 sys- 
tem and the z = 0.761 system, and the relative sizes of 
the groups at 2: = 0.716 and z = 0.758 favor the inter- 
pretation that the lens is associated with the z = 0.761 
absorption system, though the data are still inconclusive. 

4.2. Lensing Effects from the Environment 

We model the groups as SIS mass distributions and 
find the total group contribution to the convergence to 
be Kgrps ~ 0.08 (Table [4]). The addition of more group 
members could increase or decrease the group velocity 
dispersion estimate, and we therefore assume that our 
current data provide a reasonable estimate of the true 
velocity dispersion of the groups. The biased sampling of 
group members (i.e., excluding potential members close 
to the lens) further complicates our ability to assess the 
group contributions to the lensing. The lack of spectral 



coverage around the lens would tend to inflate our esti- 
mate of the group offsets from the lens; the convergences 
reported in Table [3] may therefore be underestimates be- 
cause Ksis oc 1/r. However, neither the lens nor any 
galaxies excluded from our spectroscopy are candidate 
brightest group galaxies that might be expected to lie 
near the center of the group potential; we therefore do 

not expect this bias to greatly affect our results. 

In previous work we have shown (|Auger et al.l l2007f) 
that the individual galaxies tend to provide a much larger 
contribution to the environmental convergence than the 
group halo if the lens is not the brightest group galaxy, as 
is the case with SBS1520. Therefore, all galaxies within 
a projected distance of several truncation radii from the 
lens need to be accounted for to adequately model the 
lens environment. Our HST imaging shows 5 galaxies 
within a projected distance of 80 h^^ kpc (for z = 0.761) 



Lens Galaxy Properties of SBS1520+530 
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TABLE 4 
SBS1520 Groups 



Mean Position 




Lens Offset 




(7 






RA 


Dec 


N 7 

gals 


(") 


z 


km s^-"- 


(h-i Mpc) 




15 21 48.16 


52 53 19.5 


5 


94 


0.716 


121 ± 32 


0.92 


i 0.005 


15 21 38.22 


52 55 23.1 


13 


69 


0.758 


465 ± 86 


0.45 


0.02 


15 21 43.37 


52 54 52.0 


10 


14 


0.818 


369 ± 118 


0.87 


0.06 



TABLE 5 
Group Members 



Group Redshift 


Galaxy Name 






R.A. 






Dec 






z 






0.716 


J152141. 516+525422 


87 


15 


21 


41 


516 


52 


54 


22 


870 





7154 


21.317 


21.806 




J152201. 879+525154 


94 


15 


22 


01 


879 


52 


51 


54 


943 





7159 


20.258 






J152125.755+525626 


37 


15 


21 


25 


755 


52 


56 


26 


366 





7166 


20.670 






J152154.236+525212 


74 


15 


21 


54 


236 


52 


52 


12 


737 





7166 




23.704 




J152157.414+525140 


39 


15 


21 


57 


414 


52 


51 


40 


395 





7170 




23.955 


0.758 


J152149.639+525357 


38 


15 


21 


49 


639 


52 


53 


57 


385 





7539 


21.180 


22.174 




J152128.545+525631 


31 


15 


21 


28 


545 


52 


56 


31 


306 





7539 


21.450 


22.145 




J152130.843+525642 


03 


15 


21 


30 


843 


52 


56 


42 


026 





7548 


20.740 


21.215 




J152152.978+525456 


06 


15 


21 


52 


978 


52 


54 


56 


062 





7572 


21.494 


22.418 




J152136.775+525504 


75 


15 


21 


36 


775 


52 


55 


04 


751 





7573 


21.562 


22.357 




J152133.004+525614 


24 


15 


21 


33 


004 


52 


56 


14 


241 





7575 


21.118 


22.169 




J152148.084+525229 


23 


15 


21 


48 


084 


52 


52 


29 


227 





7590 




23.202 




J152130.631+525643 


55 


15 


21 


30 


631 


52 


56 


43 


555 





7593 


20.882 


21.878 




J152140.236+525420 


07 


15 


21 


40 


236 


52 


54 


20 


069 





7594 




23.251 




J152128.163+525657 


95 


15 


21 


28 


163 


52 


56 


57 


952 





7594 


21.088 






J152133.198+525640 


07 


15 


21 


33 


198 


52 


56 


40 


068 





7599 


20.984 


22.100 




J152144.852+525448 


52 


15 


21 


44 


852 


52 


54 


48 


521 





7604 


17.995^= 


18.071<= 




J152139.952+525435 


53 


15 


21 


39 


952 


52 


54 


35 


532 





7630 




23.862 


0.818 


J152139.850+525524 


08 


15 


21 


39 


850 


52 


55 


24 


082 





8132 


22.137 


23.352 




J152205.728+525538 


52 


15 


22 


05 


728 


52 


55 


38 


522 





8168 


21.180 


21.663 




J152153.761+525317 


33 


15 


21 


53 


761 


52 


53 


17 


327 





8172 


20.504 


22.846 




J152148.911+525708 


53 


15 


21 


48 


911 


52 


57 


08 


526 





8175 


21.901 


23.341 




J152128.683+525352 


65 


15 


21 


28 


683 


52 


53 


52 


653 





8178 




23.690 




J152125.164+525411 


23 


15 


21 


25 


164 


52 


54 


11 


229 





8178 




23.096 




J152130.942+525514 


19 


15 


21 


30 


942 


52 


55 


14 


195 





8182 


21.547 


21.961 




J152146.267+525547 


69 


15 


21 


46 


267 


52 


55 


47 


690 





8185 


20.915 


21.823 




J152154.462+525233 


71 


15 


21 


54 


462 


52 


52 


33 


711 





8197 




23.733 




J152139.973+525532 


23 


15 


21 


39 


973 


52 


55 


32 


229 





8226 




23.011 



^ From SDSS imaging <Adelman-McCarthv et al.|[2007fl .'^ From J. P. Kncib, private communication.^ Magni- 
tudes include emission from the lensed quasar. 



of the lensing galaxy (Figure [3]) , and the colors of these 
galaxies are consistent with being at approximately the 
same redshift as the lens galaxy. Four of these galax- 
ies have not been treated in previous lens models, and 
assigning all four galaxies a redshift z — 0.76 and veloc- 
ity dispersion a = 120 kms~^ yields an additional con- 
vergence of Kgais ~ 0.11 for SIS profiles. Previous lens 
models fo r SBS1520 have in cluded cluster halos in the 
modeling (Bur ud et al.ll200^ . but the contribution from 
these four individual galaxies was not accounted for in 
those models and can therefore be treated as a correc- 
tion to previous models. 

4.3. Modeling the Mass Slope of SBS1520 

Early-type galaxies have been shown to have nearly 
isot hermal mass profi l es wit h only a small intrinsic scat- 
ter (|Koopmans et al.l l2006l ) and it is therefore reason- 
able to use SIS mass distributions to model lenses with a 
small num ber of constrai nts. However, part of the scat- 
ter in the iKoopmans et al.i , (2006) result may be due 
to systematic effects (jAuger 2007) and our improved 
knowledge of the lens environment and the ellipticity of 



the lens galaxy makes it worthwhile to create updated 
lens models for the system that do not re quire the as- 
sumption of an isothermal mass profile. iDobke et al.l 
(2007) have used N-body simulations to suggest that 
galaxies in overdense environments may undergo strip- 
ping that creates steepened mass slopes when interact- 
ing with other galaxies. This presents a physical mo- 
tivation for allowing for non-isothermal profiles in lens 
modeling and may explain why some len ses have system- 
atically underestirnated i Jn (e.g.. (Kooprnaris et al.|[2000t 
iTreu fc Koopmansl[200l lAuger et al.ll2007t l. Consider- 
ing that several independent estimates of Hq have con- 
verged to approximately the same value of approximately 
72 km s~^ Mpc~^, we have chosen to fix the value of Hq 
at 72 km s~^ Mpc~^and allow the mass slope to vary in 
our model fitting. 

We model the galaxy groups as SIS halos with the 
velocity dispersions and positions reported in Table [4] 
and we model all galaxies within 80 h~^ kpc (ex- 
cept Galaxy M) as SIS halos with velocity dispersion 
a = 120 kms~^. We model Galaxy M as an SIS halo but 
leave the velocity dispersion as a free parameter. The 
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Fig. 7. — The spatial distribution of galaxies in the groups at (a) z = 0.716, (b) z = 0.758, and (c) z = 0.818. The location of the lens is 
marked with a dashed crosshair and the circle is centered on the group position and has a radius of 1 h~^ Mpc at the redshift of the group. 



lens is modeled with lensmodel (|KeetoDll200lD as a power- 
law mass distribution, p oc r~", and the slope and lens 
strength are allowed to vary in the optimization; the el- 
lipticity, position angle, and location of the lens are fixed 
to the values obtained from our AO imaging. T he time 
delay is At = 130 ± 3 days (|Burud et alj|2002[) and, in 
principle, the flux ratio is best determined by using the 
time-delay co rrected ratio from monitorin g campaigns 
A/B = 1.89. lBurud et al.ll2002i) . However. iBurud et all 
2002f ) do not perform their modeling with this value and 
do not directly comment on why the time-delay flux ratio 
was ignored; we were also unable to produce models of 
the lens that converged with such a low flux ratio. We 
use the quasar Mgll emission line from our LRIS data to 
determine the flux ratio to be 2.9±0.4, where the error is 
derived from the scatter between the seven independent 
LRIS spectra. We choose the quasar Mgll line because it 
should be relatively insensitive to quasar variability and 
it is the reddest line available, th us minimizing proble ms 
with differential extinction (e.g.. iWisotzki et al.ll2003l ). 

SBS1520 can be reasonably well modeled with a 
steeper than isothermal power law profile {a = 2.29^0'"^) 
if -ffo is set fixed to 72 km s'^ Mpc'^ (Table O; with 
three degrees of freedom (DOF), xVDOF = 1.1. The 
main contribution to the chi-square comes from the flux 
ratio; the lens favors a flux ratio higher than the ratio 
derived from our data. We note that Burud et al. (200^ 
used a flux ratio of 3.6 derived from emission lines in 
their spectra, though this is substantially higher than 
the emission line ratio from our spectroscopic data and 
the continuum flux ratio of 2.4 that we find in our NIRC2 



K-band imaging. 

There is a well-kn own degeneracy between the mass 
slope and Hq (e.g., IWiUiamslF Sahal l2000t iWucknitd 
l2002f ) that allows a substantial amount of play between 
these parameters (compare our results with IBurud et al.1 
120021 for example). The value of Hq obtained for 
SBS1520 for an isothermal mass distribution when the 
environment is fully taken into account may be as low as 
Hq k 46 km s~^ Mpc~^if we account for the unmodeled 
convergence due to the environment, Kext, via the simple 
relation 

Ho true ^ (1 ^ext)HQ „iodel 

and use the isothermal models of Buru d et al.l (|2002f ). 
The mass slope modeled in our analys is is steeper than 
expec ted for early-type galaxies (e.g., iKoopmans et al.1 
l2006f ). but the presence of a nearby galaxy, Galaxy M, 
suggests that SBS1520 may be undergoing interaction- 
induced steepening fe.g.. lDobke et al.ll2007l ). iRead et all 
|2007) have used pixelated mass models to model 
SBS1520 and find a nearly isothermal profile, a = 1.95. 
These models are able to account for 'shape degenera- 
cies' that may also ef f ect in terpretations of time de- 
lays dS aha fc Williams' ^ 20061 ). though the analysis of 
iSaha fc Williams (2006) indicates that shape degenera- 
cies are probably not important for SBS1520. Note that 
iRead et al] ()2007[ ) do not explicitly model the lens envi- 
ronment and they are therefore measuring the joint pro- 
jected mass slope of the lens and environment. If the 
lens is not near the (projected) group center, th inferred 
mass slope would be shallower than the true mass slope. 
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TABLE 6 
Lens Model Parameters 



Parameter 



Value 



Lens Position 


-1('1531, 


-0^4150 


QSO A Position 


0"0, O'.'O 




QSO B Position 


-1'.'4227, 


-0'.'6569 


Flux Ratio (A/B) 


2.9 




At 


130 days 




Ho 


72 km s" 


-1 Mpc-i 


Power Law Index'' 


2.29 




Lens Strength'' 


0.53 




f GalaxyM^^ 


84km s^ 


1 




3.4 





^ Parameter was varied in the model fitting. 



5. CONCLUSIONS 

We have obtained deep AO imaging and optical spec- 
troscopy of the time-delay lens SBS1520. The AO imag- 
ing has allowed us to fix the lens galaxy properties with a 
high degree of precision when performing the lens mod- 
eling, and the data indicate that the lens has an elliptical 
morphology and perhaps a disk. The new spectroscopic 
data suggest that previous determinations of the lens red- 
shift may be incorrect, and the data also allow us to 
quantify the lensing contribution of several groups in the 
immediate foreground and background of the lens. Lens 
models created with these new data can be well-fit with 
a steeper than isothermal mass slope (a = 2.29) if Hq 
is fixed at 72 km s^^ Mp c~^; isothermal mo dels require 
Ho ~ 50 km s^i Mpc^^. iDobke et all ()2007f ) found that 
galaxies in overdense environments might have steeper 
than isothermal mass sl opes caused b y interactions with 
other galaxies (also see lAuger1l2007D . This suggests an 
interpretation that we are observing transient steepen- 
ing of the mass profile due to galaxy-galaxy interactions 
and indicates that other lens systems that have obtained 
anomalously low values of Hq may lie in overdense re- 
gions and near an interacting galaxy (e .g., B1600-I-434; 
iKoopmans et al.ll2000l : lAuger et al.ll2007f) . Alternatively, 
SBS1520 can be modeled in a manner consistent with 
an isothermal profile and Hq = 64 km s'-^ Mpc~^ if 
if the lens is modeled by a pixelated mass distribution 
and .j ointly modeled with other lens systems (|Read et al.l 
|2007( ). These models indicate that twisting ellipticity, 
triaxial projection effects, or other shape degeneracies 
may be ef fecting the parametric analyses of SBS1520 
(|Saha fc WiUiam s 2006). 

However, there are still several ambiguities in the data 
that need to be resolved before making definitive claims 
about the profile of SBS1520, particularly in the context 
of the interaction-induced steepening scenario. While we 
have argued that the lens redshift is likely to be z = 0.761 
and not z = 0.717, the data are not conclusive. Further- 
more, our modeling has assumed that all of the neighbor 
galaxies are at the group redshift; if this is the case, the 
z = 0.76 group centroid would be pulled closer to the 
lens and the group would therefore provide a larger con- 
tribution to the lens model. If this is not the case, the 
neighboring galaxies might have a smaller impact on the 
lens model. This is particularly important for Galaxy 
M, as this is the neighboring galaxy that most affects 
the lens model but also has colors least like the lens 
galaxy compared to the other field galaxies. It is also 



important to verify that at least one of the neighboring 
galaxies is at the same redshift as the lens because this 
is a requirement of the interaction-driven steepening hy- 
pothesis. Finally, obtaining a dynamical estimate of the 
lens mass would help to further constrain models and 
potentially distinguish between shape degeneracies and 
the mass slope degeneracy. 
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